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Ch . Abstract 

1—j | The binding of a system of ./V polarons subject to a constant magnetic field of strength 

jy-j ' B is investigated within the Pekar-Tomasevich approximation. In this approximation the 

O^l ■ energy of N polarons is described in terms of a non-quadratic functional with a quartic 

term that accounts for the electron-electron self-interaction mediated by phonons. The size 
of a coupling constant, denoted by a, in front of the quartic is determined by the electronic 
properties of the crystal under consideration, but in any case it is constrained by < a < 1. 
For all values of N and B we find an interval ajv,s < ot < 1 where the N polarons bind in a 
single cluster described by a minimizer of the Pekar-Tomasevich functional. This minimizcr 
is exponentially localized in the ./V-particle configuration space M. 3N . 

<N ! 1 Introduction 

> ■ 

The electron-phonon interaction in a polar crystal mediates an interaction between pairs 
\Q \ of electrons which becomes an electrostatic Coulomb attraction in the Pekar-Tomasevich 

■ approximation. This attraction competes with the Coulomb repulsion between the equally 
charged electrons, and the question arises whether TV electrons may form a bound cluster. 
Due to the constraint on the parameters of the model, the l/|x|-part of the electron-electron 
interaction is repulsive. There remains, however an attractive short range interaction, which 
seems to be of van der Waals type and which may lead to N— particle bound states [11] . 

• i-h ' This phenomenon of bound multipolarons had previously been observed in Frohlich's large 

■ polaron model on which the Pekar-Tomasevich approximation is based [171 0] . Similarly, 
, the binding of polarons subject to a constant magnetic field had been investigated within 

the Frohlich model [3! . Yet, in that case, the analysis in the literature is based on poorly 
justified variational estimates, and the conclusions remain doubtful. The present paper 
establishes, within the Pekar-Tomasevich approximation, the existence of bound iV-polaron 
clusters in a constant magnetic field of any strength. It is a continuation of a previous work 
of one of us, concerning the case N = 2 [5]. 

The Pekar-Tomasevich approximation to the large polaron model of Frohlich describes 
the energy of N polarons through an effective functional that depends on the wave function 
^ € Wn '■= A N L 2 (R 3 x {1, . . . , q}) of the particles only. We are mainly interested in the 
case of spin- 1/2 fermions but we can allow for arbitrary q £ N without more effort. The 
functional is then given by 

JV 
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where U, a > are constants, and 

N q . 

P®( x ) : = X) X / l*fel'---'£j-l'( a; ' Cr j)'£j + l'---'2jv)| 2rf £l---^j--- d 2Ar, (2) 

is the density associated to VF We have introduced the notations Xj = (xj,crj) for elements 
of R 3 x {1, . . . , q} and we set J dx_j = Y^j =i I ^ x r Of course in ([2]) the sum with respect 
to j may be replaced by a factor of N, due to the symmetry of but we shall allow for 
Boltzons later on, and hence we prefer (J2j as the definition of p^. Furthermore, Da,x '■— 
— iW + A(x) where the vector potential A : M 3 — > M 3 generates a magnetic field B = curl A 
We are primarily interested in the case where B is constant and hence A will be assumed 
linear. The positive parameters U, a are constrained by a < U due to their role in the 
Frohlich large polaron model. Mathematically, any real values are conceivable for U and 
a, but < U < a leads to thermodynamic instability [10] . The energy of the fields 
U 3N/2 ^(Ux 1 ,cr 1 ,...,Ux N ,a N ) and UA(Ux) upon the substitutions Ux — > x and a/U —> a 
becomes proportional to U 2 . We therefore set U = 1 and we require that < a < 1. 

It is easy to see, using the diamagnetic and the Hardy inequalities, that £ N ' a is bounded 
below if restricted to the unit sphere ||\&|| = 1. The minimal energy, 

:= inf £ Na (^) 1 (3) 
F1 ll*ll=i 

is therefore finite. By moving particles apart, one can see that Ep^ < E P \^ + Ep T k ' a for 
k = 1, . . . , N — 1. The question is, whether it takes energy to do this, that is, whether for 
some a < 1, 

i [E k p % + E^ T k ' a } - Eprp~* > 0. (4) 
Our main result is the following theorem: 

Theorem 1.1. Assume that the vector potential A is linear (constant magnetic field B). 
Then, for all N G N there exists un^b < 1 such that for aj^^s < a < 1 and U = 1: 

(a) the binding inequality (j4|) holds, 

(b) the functional {1} has a minimizer. 

Analog results hold in the case of bosons and boltzons, that is, for Hn = L 2 (M. 3 x 
{1, . . . , q}), the symmetric product of N copies of L 2 (R 3 x {1, . . . , q}), or H N = <E) N L 2 (R 3 ) 
without symmetry requirements. The proofs in these cases are similar and in the case of 
Boltzons the proof of (a) becomes much easier. Yet the property (a) even for boltzons is 
a subtle correlation effect since the restriction a < 1 means that the Coulomb repulsion 
dominates the attraction for states of the form (p\ ® . . .® tp N . We remark that Theorem ll.il 
has consequences for the binding of boltzonic polarons in the large polaron model of Frohlich 

Hi- 

For a = there is no minimizer and, in the absence of magnetic fields, there is no 
binding for a small enough 6j. The existence of a minimizer is a phenomenon due to the 
non-linearity and it occurs whenever the binding inequality (j4|) is satisfied (and a > 0). 
For other non-quadratic energy-functionals associated with many-body quantum systems 
this has previously been pointed out and described as a non-linear HVZ-Theorem [111 [7] . 
In this paper we show that (a) (6) is a consequence of a linear HVZ-Theorem for an 
iV-body Hamiltonian that is intimately related with the physics of the polaron problem: 
there is a Hamiltonian H a depending on a charge density a € i 1 (R 3 ) such that 



/XEprp := mm 

rJ Kk<N- 
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with equality for a — pqj. We may think of aa as the charge density due to a hypothetical, 
possibly non-optimal, lattice deformation caused by the electrons. For £ N ' a (^ n ) near Ep T , 
(\l/ n ) being a minimizing sequence with densities (p n ), the binding inequality implies that 
H Pn has an isolated ground state energy separated from the essential spectrum of H Pn by a 
gap that is uniform in n along a subsequence. This uniformity implies uniform localization 
of (or concentration of minimizing sequences) up to magnetic translations. 

Our proof of part (a) in Theorem II .11 is based on a variational argument that is inspired 
by [8] but is considerably more involved in the present case of particles with statistics. 

The following theorem gives further information about the minimizers found in Theo- 
rem [TTT] In Theorem 11.21 and throughout the paper we use the notation V p := p * | • 

Theorem 1.2. Iffy 6 T~Ln *s a minimizer of£ N,a , then it solves the non-linear Schrddinger 
equation 

1 \ 

k=l i<j 1 1 3 1 J 

where A € K is the lowest point in the spectrum of the Schrddinger operator on the left 
hand side and p is the density of ^ . Moreover, if (j4|) holds then the spectrum of the 
Schrddinger operator on the left hand side is discrete below A + AEp^* and hence if (3 6 M 
with P 2 < AEp^, then 

e^ 1 * £ H N . (6) 

In the case N — 1, A = the Pekar-Tomasevich functional reduces to the Pekar or 
Choquard functional which is well known to be minimized by a spherically symmetric, 
positive function that is unique up to translations [T2l US] . 

Existence of a magnetic polaron and the binding of two polarons subject to an external 
magnetic field, not necessarily constant, was previously established in 8 . In the present 
paper, the methods developed in [8] are extended and generalized to the case of N > 2 
particles of fermionic, bosonic or bolzonic nature. Results similar to ours in the case A = 
where previously obtained by Lewin in [TTJ. Lewin establishes a bound on the binding 
energy of the form of a van der Waals potential with exponentially small corrections. To 
this end he uses the variational state introduced by Lieb and Thirring in connection with 
the van der Waals binding of neutral atoms and molecules |15j . This approach makes crucial 
use of spherical averaging and the Newton's theorem. It brakes down in the presence of a 
magnetic field where the rotational invariance of £ N < a is broken. Moreover, in the absence 
of a magnetic field our Theorem 11.21 gives more information than the corresponding result 
of Lewin, as it relates the binding energy AEp^ to the gap between A and the essential 
spectrum of the Hamiltonian in ([5]). Lewin, in the case of binding, merely finds that such 
a gap exits provided that a > 1 — 1/N. 

The Theorem 11.21 opens the following new view upon the phenomenon of A^-polaron 
binding: if a Hamiltonian of the type in ([S]) with total positive charge aN is shown not to 
bind N electrons, then binding of N polarons is excluded. Here binding means positivity 
of the binding energy. - In the case where the density is spherically symmetric and A = 
we deduce from [13] that the Hamiltonian has no ground state if a < (1 — N~ 1 )/2. This 
leads to the following corollary: if a < (1 — N~ 1 )/2 then a hypothetical minimizer of £ ' a 
cannot have a spherically symmetric density. 

This paper is organized as follows: In Section [5] we outline the proof of our main 
Theorem and we introduce the most important tools. In Section [3] we prove an operator 
inequality which is of crucial importance for the proof of existence of a minimizer of the 
Pekar-Tomasevich functional, as well as the proof of the second part of Theorem 11.21 In 
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Section |4] we use the operator inequality to prove existence of a minimizer and exponential 
decay of any minimizer of the Pekar-Tomasevich functional. In Section [5] we establish the 
binding inequality 

Acknowledgements. The first author (I. A.) is grateful to Fabian Hantsch and David Wellig 
for numerous stimulating discussions and for introducing him to the theory of multipo- 
larons. He also thanks Mathieu Lewin for interesting discussions on iV-body quantum 
systems. 

2 Preparations and elements of the proofs 

The minimal energy Ep£ is continuous in a because it is concave in a as the infimum of the 
afiine functions a h- > £ k > a (^>). Hence, it suffices to establish the binding in the case a = 1. 
Our proof that binding implies existence of a minimizer, i.e (a) ==>■ (b) in Theorem 1 
as well as the proof of Theorem 11.21 readily generalize from the case a = 1 to any a > 0. 
We therefore put a = 1 for notational simplicity, that is, 

n*>:=(^f>L,+g^)*)-^>. (7) 
where D(p) := D(p,p), 

D M :=U<P^ dxd y, (8) 
2 J \x-y\ 

and 

E":= inf £ N (V). (9) 
ll*ll=i 

The domain of £ N is the form domain, Qn,a, of Ylk=i ^A x k -> tnat Qn,a = G Hn ■ 
DA,x k ^ € L 2 ,Vfc <E {1, ...,N}}, and we use || • \\q N:A for the corresponding form norm. By 
a minimizer of £ N we shall always mean a normalized vector ^ € Hn with "J € Qn.a and 
£ N (^) = Ep T . Throughout the paper we use (•, •) and |j • || for the usual inner products 
and norms of (g) 7V L 2 (R 3 x {1, . . . , q}) and Hn- 

By the above explanations it remains to prove the following theorem in order to establish 
Theorem 11.11 

Theorem 2.1. Assume that the vector potential A is linear. Then, 

(a) there exists a minimizer of £ l , 

(b) if £ 1 J £ N_1 have minimizers then 

E^ T <E PT + E^ k , Vfc = l,...,iV-l, (10) 

(c) if (|10j) holds then £ N has a minimizer. 

Part (a) of Theorem l2.1l is known from [8] but we shall reprove it as a part of the proof 
of part (c). Part (b) is proved in Section [5] by variational arguments. Sections [3] and |4] are 
devoted to the proof of (c). The remainder of the present section describes the difficulties 
met in the proof of (c) and collects our tools for dealing with them. 

Any proof of (c) must deal with the following translation invariance of £ N : Let A : 
R 3 M 3 be linear, h G M 3 , and * e H N . If T h ^> is defined by 

JV 

(T h ty)(xi, ai, . . . , x N , ctjv) = Y[ e'^^^ix! +h,ai,...,x N + h, a N ), (11) 

3=1 
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then pT h w(x) — p^{x + h) and 

e N (T h *)=£ N (9). (12) 

Due to (llip and (|12[) a minimizing sequence of £ N may converge to the zero function 
weakly. On the other hand in view of Lemma \A.1\ a weak limit ^ £ Hn with = 1 is, 
indeed, a minimizer of £ N . Our task is thus to find a minimizing sequence of £ N that does 
not suffer any loss of norm in the limit. One of our tools to this end is the following form 
of the Concentration Compactness Principle [TB] : 

Proposition 2.2. Let (pk)k>l be a sequence of nonnegative functions in L 1 (M 3 ) with 
J pk = N. Then there exists a subsequence of (pk), denoted by (pk) as well, such that one 
of the following holds: 

(i) (Vanishing) For all R > we have that 



lim sup / pk = 0. 

k ^°°yem 3 Jb( v ,r) 

(ii) (Dichotomy or compactness) There exists A £ (0, N] such that for all e > there exist 
R e > 0, a sequence yu = yk(s) in R 3 , and a sequence P& = Pfc(e) in R with Pk — > oo 
as k — > oo, such that the sequences of functions^ 

Pk.l ■= PkXB(y k (e),R e ) 
Pk,2 ■= PkXB(y k (e),P k (e))C 

satisfy for k > ko(s) the bounds 

\\pk - pk,i - Pfe,2||z,i < e, (13) 
|||pfc,i|Ui-A| <e, |||p fc , 2 |Ux-(iV-A)|< £ (14) 

and 

dist(s\ipp Pfc,ij supp pk,i) — ^ oo, {k ->• oo). (15) 

If m is a positive integer such that mX > N, then after passing to a subsequence once 
more, there exists s±, ...,e m -x > 0, and S > such that 

liminf / pk,i > 5 (16) 

/or all e > small enough. Here R>k,s — B(yk(s), R £ ). 

Proof. We shall only prove the last part of (ii) . The rest is a variation of the Concentration 
Compactness Principle. Let Ex, Si > be such that m(A — ei) > N + mSi. Assuming that 
the lemma is wrong we inductively construct E\ > £2 > ... > s m > and a subsequence of 
Pk denoted by pk as well, such that 



Pk < J—J' ^ l = 2 ' •■•' TO - 



Using this together with (TT4"]) and the inequality 



7 — 1 2<j<771 



denotes the characteristic function of the set A and B(y, R) is the ball of radius R centered at y in '. 
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we obtain that 



~ m 

liminf / pk > T^(A - £j) — (m — l)6i > m(\ - ei) - (m — l)Si > N, 

where the last inequality follows by the choice of S\ and 8\. This is in contradiction with 
/ Pk = N, which concludes the proof of the lemma. □ 



The following lemma is the reason for the new part (|16j) in the above version of the 
Concentration Compactness Principle. 

Lemma 2.3. In the case (ii) of Provosition UTM if (pk) * s chosen to satisfy (|16p, then 

liminf D(p K1 ) > (17) 

uniformly for small enough e. (Recall that p kt i depends on e.) 

Proof. Let e be small enough for (fTB|) and fixed. By (TT5J) there exists k (e) such that for 
all k > ko(e) we have 

f (m-l)S , , 

/ Pk.i > ± '-■ (18) 

This means that f B ^ pk,i > 6/m for some j £ {1, ...,m — 1} depending on k > fco(e). 
Since diam(i? / t. ei ) = 2R £j , we conclude that 

D(p k<1 ) > I j P kAx)dx\ > miE-^, (19) 

which proves the lemma. □ 

We want to construct a minimizing sequence (^ k ) that is concentrated near the origin 
(after translations). Applying the Concentration Compactness Principle to l^fcl 2 would not 
work, because the Pekar-Tomesevich functional is invariant under translations of the form 
(jlip . only, and not under general translations in M. 3N . Thus, we apply the Concentration 
Compactness Principle to the densities, where dichotomy may mean various things for 
the wave function. Rather than trying to exclude all of them we show directly that non- 
vanishing of the sequence pk, leads to concentration of a subsequence of ^> k . This is possible 
thanks to an HVZ-type operator inequality for the Hamiltonians H^ k defined as follows: 
for a given real- valued density a G L^R 3 ) n L 6 / 5 (R 3 ) we define 

V CT := a * i (20) 

and 

N 1 
H? := " V a (x 3 )) + £ -— — + D(a), (21) 

j=l i<j 

which is well defined by the choice of a ([2] Corollary 5.10). In all the following this 
operator is considered defined in Hn unless explicitly stated otherwise. The following 
lemma, taken from [6], relates the Pekar-Tomasevich functional to the linear Hamiltonian 
USD: 
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Lemma 2.4 (Linearization of the Pekar-Tomasevich functional). For any density 

< <#,Ff#) (22) 
equality if and only if a = p<$, . In particular, for all N € N, 

ijf > Eprp. (23) 

In particular, if (^? k ) is a minimizing sequence for £ N and (p k ) is the sequence of the 
corresponding densities, then 

]im(^ k ,H^k)=E^ T . (24) 

Proof. By the definitions of ,£ N , V a and D we have that 

(% H?V) - £ N m = D(a) + D(p 9 ) - 2D(p^,a) = D(a - pv) > 0, 

where the last inequality follows from the positivity of the Fourier transform of |.| . This 
proves (JUl). Inequality (j2"3")) follows from (|2"2")l and from the definition, Equation ©, of 
E$ T . Equation ([24]) follows from (f k ,H^ k ) = £ N (^ k ) and from the choice of (* fc ). □ 

The main steps in our proof of part (c) of Theorem 12. II are as follows: 

Step 1 is to exclude vanishing for the sequence of the densities (pk) associated with a 
minimizing sequence To this end we prove that vanishing implies that D(p k ) — > 

which is easily seen to be in contradiction with £ N (^ k ) — > Ep T . 

As vanishing has now been excluded, the second alternative of Proposition 12.21 must 
apply to the densities (p k ) of any minimizing sequence (^ k ). Upon the translations fy k — > 
Ty k ^ k , see (jlip . we may assume that some part of the densities p k is concentrated near 
the origin. 

Step 2 is the proof of the operator inequality 

ff£>Egr + d(l-J e ) + 0(V£), (25) 

where d > 0, J e is compactly supported and < J £ < 1. The proof of (|25|) is based on 
the properties of p k as described by Proposition 2.2 (ii), on Lemma [2.31 and on a suitable 
partition of unity that is adjusted to the supports of p k> x and p kt 2- 

Step 3 is to show that (|2~5|) implies concentration of (^ k ). This is easily done with the 
help of ([M]) and the fact that e in (|2~5|) may be taken arbitrarily small. 



3 Absence of vanishing and the operator inequality 

Our goal in this Section is to establish absence of vanishing of the sequence of the densities 
(p k ) associated with a minimizing sequence (^ k ) and to prove the operator inequality of 
Proposition 13.21 

Lemma 3.1 (Absence of vanishing). The sequence of the densities (p k ) associated with 
a minimizing sequence (^ k ) of £ N cannot be vanishing. 

Proof. We shall derive a contradiction from the assumptions that (^ k ) is minimizing and 
that (p k ) is vanishing at the same time. The vanishing of (p k ) implies that 

lim D(p k ) = 0, (26) 

k—yoo 
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as we will prove shortly. By ([7]) and (|26[) wc have that 



lim S N (%k) > Nirio-(Dl) > N\B\. (27) 

k— >OG 

On the other hand Ep T < NEp T by general principles and Ep T < \B\, by [5]. It follows 
that 

E^ T <N\B\, (28) 

which we combine with (|27|l to conclude that the sequence (^ k ) is not minimizing in 
contradiction to our assumption. 

We now turn to the proof of (f2T)|) . From ||/?fc||Li = N it follows that, for any r > 0, 

D( Pk )<{ eMp^l dxdy+ ^. (29) 

7|x-*i<r \ x -y\ r 



and 



"^M-daafo < N sup / (30) 

\x — y\ <r |** 1 



For each x E K , by Cauchy-Schwarz, 

T^W/ PHtody) 1/2 (f T^dyY 2 . (31) 

|x-3/|<r F _ 2/1 \J|x-y|<-r / \-/|x-!/|<r \ X ~V\ J 

On the right hand side of (I31[) , the first factor vanishes uniformly in x in the limit fc — > oo, 
by the assumption that (p&) is vanishing. The second factor is bounded uniformly in x 
because of Lemma IA.1I and the estimate 

f Pk(y) , f l*fefe,---,2jv)| 2 . 7 ^,||,T, 112 

J |^p^ = Ly ^i--.^<4||vM Q ^. (32) 



Here we used the Hardy and diamagnetic inequalities. As we have now shown that ()31|) van- 
ishes uniformly in x in the limit k — > oo, we conclude, combining (|29[) -(|3T j) . that D(pk) — > 
as k — > oo because r > may be chosen arbitrarily large in (|29p . □ 

Proposition 3.2. Suppose that (| 10[) /io/ds and Zei (^fc) be a minimizing sequence whose 
densities pk = p* fc have the properties of Proposition [2J% (ii). Then there exists a sub- 
sequence of pk, denoted by pk as well, and a positive number d > such that for all 
e > small enough there exists a function J e € Cq°(M. 3N ; [0, 1]) symmetric with respect to 
exchange of particle coordinates, such that for all k > ko(e) 

H» > E$r + d(l - r yfc J e ) - N{2^C + eiV) - 2 N (eN) 2 , (33) 

w/iere y fc = y k {e) is given by Proposition\2J& ( ii), T yk J e (xi, x N ) := J s (xi—yk, — , x N -y k ) 
and C :— 2 sup ||^fe|| g N A < oo ('see Lemma \A.l\) . If the sequence (pk) is concentrated, i.e. 
if A = N in Provosition \2. S\ (ii), then we may choose d — AE N := min{_Ep T + Ep^ k \ k — 
1 V !|- E^ T . 

We fix e > and (^k) as described in Proposition 13.21 Let (yk) be the corresponding 
sequence provided by Proposition 12.21 (ii). After the translations ^k H> T yk ^k defined by 
Equation (fTTj) we may assume that the densities of (^k) have the properties of Proposition 
12.21 (ii) with yk — 0. It thus remains to prove Proposition 13.21 in the case yk = 0- As a 
preparation we will first establish the following two lemmas. 
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Lemma 3.3 (Partition of unity). Let e and ^k be as explained above. Let also pk = p>s k 
and pk.i be as in Proposition \2.2\ (ii). Then there exist kg > 1 and non-negative functions 
jxj 2 ■ R 3 -> R with 

0<ji,j 2 <l, J?+jI = l, ||Vj*||i- < e, mippiiCB(0,Ji e + |), (34) 



such that for all k > ko, 



dist(snpppk,i,swppj 3 -i) >-, i = 1,2. (35) 

e 



If a = (ai,...,ajv) € {1,2} , i/ien the functions 

N 



J a (xt,...,XN) ■= Y[j aj (xj) (36) 



/icwe i/ie following properties: 



0<J„<1, J! J « = 1 ' HVJalli- <£^V. (37) 

aG{l,2}« 

Proof. It is an elementary exercise to construct non- negative functions /i , /2 € C 100 (R) 
with sup.,, 1/^ (x)| < 1, ff + /| = 1, /i = 1 on (—oo, 1] and / 2 = 1 on [3, oo). Let 

Je (x) = M(\x\-R £ )e). 

Using the properties of fi, f 2 and the fact that -Pfc(e) > R e + 4e~ 1 for k large enough, see 
Proposition 12.21 (ii), one easily verifies that ji,j 2 have the desired properties. (|3"?| follows 
from (|36[) and the properties of j\ , j 2 ■ □ 

Lemma 3.4. Le£ e and (^k) be as in Lemma \3.3l and C :— 2sup ||^fc||Qjv A as in Propo- 
sition [3~M If Pk,Pk,i are given by Provosition \2.2\ (ii), then for k large enough, 

V Pk ~V PkA ~V PK2 < V~eC, (38) 
(V Ph -V Pk<i )jf < (V^C + eN)jf, i = 1,2. (39) 

Proof. By the definitions of V Pk , V Pk 1 , and V Pk 2 , we have 



V Pk - V Pk ,i - V Pk ,2 = {Pk ~ Pk.l - Pk,2) 



* 



where < pk — Pk.i — Pk,2 < Pk- Hence, by Cauchy-Schwarz, (f3"2j). and ([To 

\(V Pk - V Phtl V Pk , 2 ){x)\ < (| J^^dy) ' (/ iPk - Pk,i - Pkfl)dyj ' 

To prove (|39 )l . by (|38 )l it suffices to show that V Pk 3 _jf < eNjf. This easily follows 
from Qg) and \\p k ,3-i\\v < \\pkWv = N. □ 

Proof of Proposition \3.2l In this proof we shall tacitly assume that k is large enough 
so that the statements of the previous lemmas apply. By the IMS localization formula [5] , 

K= E JaBS,j a - J2 i VJ "i 2 - ( 4 °) 

oe{i,2}« oe{i,2}™ 
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We will now estimate the terms J a H^ h J a from below. 

1st Case : a has n ones and N — n twos, < n < N. We may assume without loss of 
generality that a = (1, 1, 2, 2). From p k > Pk.i + Pk,i it follows that 

D( Pk ) > D{p ka ) + D(p k , 2 ). (41) 

This, together with (131)1) and (|2"5|) implies that 

JaH^Ja > Ja{H n pkl + H^)J a - N(^C + eN)J 2 , 

> {E n PT + E^ T n )Jl - N(V^C + eN)J*. (42) 

Note that H™ acts on the coordinates labeled by 1, .... n, while H^~ n acts on the ones 
labeled by n + 1, . . . , N, Moreover, Ja'H.N C H„ ® %N—n by construction of J a . 
2nd case : a = (2, . . . , 2), i.e., only twos. By (j39j) and (|4i"j) . 

J a H%J a > J a (D( PkA ) + H^ 2 )J a - N(^C + eN)Jl (43) 

By (|2"B")) we have ff^ 2 > -Ep T and by Lemma 12.31 there exits a constant 7 > such that 

D(pk,i) > 7, for e small enough. (44) 

It follows that, for e small enough, 

JaH N p J a > {Eprp + 7 ) J Q 2 - JV(v/iC + sN)Ja- (45) 

3rd case : a = a := (1, . . . , 1). Since > Ep T , we have 

J ao H™ k J ao > E^ T J 2 Q . (46) 

Combining the results (|4"2")l , (|4"5j) and (|4"6")l from the three cases above with (|3"T)) and (|^U|) 
we obtain (l33l) with J e = J 2 o and d = min{ 7 , AE N }, which is positive due to the binding 
assumption ([TO)) . 

In the case A = N we may improve our bound in the second case to get d — AE . 
Indeed 

N 

K,2 > E { D %*i - V p^ ^ NEp T - NC^~e, (47) 

because D 2 Ax . > Ep T and V Pk 2 (x) < CWpk^W^i < C^fe by the Cauchy-Schwarz, Hardy 
inequalities 

we conclude that 



and diamagnetic inequalities. Here we used A = N and (JUJ). Since NEp T > Ep T + Ep T 1 



JaH^Ja > (E» T + AE N )J 2 a - N(V^2C + eN)J 2 a , 



which we use in place of (|45|) . □ 



4 Existence of a minimizer and exponential decay 

In this Section we prove parts (a),(c) of Theorem 12. II and then we prove Theorem ll.2l The 
part (b) of Theorem 12.11 will be proved in the next Section. 

Lemma 4.1. Assume that (|10p holds. Then, there exists a minimizing sequence (3>fc) with 
the following property: for every 5 > there exists P > such that 

liminf / |$ fe | 2 >1-S. (48) 
fc ->°° Jb(o,p) 
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Proof. Without loss of generality we may assume that S < 1/2. By Lemma 13.11 there 
exists a minimizing sequence (^k) for which the sequence (pk) of the associated densities 
satisfies the properties of Proposition 2.2 (ii) and hence Proposition 13.21 applies to (4^). 
The operator inequality (I33[) implies that 

H%* k ) >E^ T + d- d(* k ,T Vk J e * k ) - N{2^C + eN) - 2 N (eN) 2 . 

Upon rearranging this inequality, it follows from (pM)) that 

N 2 N 
liminf(* fc ,T y)e J e * fc ) > I - —{2C^ + eN) - — {eNf > 1-5, 
fc— >oo a a 

for e small enough. Since J e is compactly supported and < J e < 1 it follows that 

liminf f |* fe | 2 >1-S, (49) 

k^oo J B (y k ,R) 

where R and y k depend on e and hence on 6. Using an argument of Lions (see [16] ) we 
shall now replace (yk) by an other sequence (y' k ) that is independent of 6 such that (|4*^|) 
still holds after enlarging R. Let i?' and (y£,) be determined in the same way as R and (yk) 
in the case 5 — 1/2. That is, 

liminf f |* fc | 2 > 1. 

JB(y' h ,R>) * 

Since H^fcH = 1 and since 1 — <5 > 1/2, by assumption, the balls B(yk,R) and B(y' k ,R') 
must overlap for fc large enough. It follows that 

liminf / |^ fc | 2 > 1- 5. (50) 

fc ^°° JB(y' k M'+2R) 

The sequence $fe = T y iJS>k is minimizing and it satisfies (|48l) with P = R' + 2R. □ 

Proof of Theorem \2.1\ (a), (c) (existence of a minimizer). Let ($&) be given by 
Lemma [4.11 By Lemma fA. 11 part (b), (&k) is bounded in Qn,a and hence, after passing 
to a subsequence, we may assume that — >• $ G Qat,a weakly in Qjv,a- Since A is locally 
bounded it follows that — > $ locally in %jv and weakly in Hat. Hence, by Lemma 14.11 
for every 6 > there exists P > such that 

1 = lim ||$a,.|| 2 > ||$|| 2 > / \$\ 2 dx = liminf / |$ fc | 2 da; > 1 - (5. 
fc^oo is(o,P) fe^ 00 Jb{o,p) 

It follows that || $|| = 1 and hence that $fe -> $ strongly in H N . Since $fc -> $ € Qjv,a 
weakly in Qat,^, the parts (a) and (c) of Theorem 12.11 follow from Lemma fA. 11 (c). □ 

Proof of Theorem This proof is based on Lemma 12.41 which clearly holds for any 
a > 0. Let ^ be a minimizer with density p. By Lemma [2~4l Hp > Ep T and ("J, H^^f) = 
Ep T . It follows that ^ belongs to the domain of the Friedrichs' extension of Hp and that 
= Ep 1 ^. This equation agrees with the Schrddinger equation ([5]) upon subtracting 
D(p)^> from both sides. 

By pQ eigenvalues of Hp below 

S := lim ( inf ($, H"<f>)) , 

R^oo \ $e£) J? , ||*|| = 1 1 ) 

D R := {$ g Q NA | = for |z| < R}, 
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are associated with exponentially decaying eigenfunctions. This means that e^ 11 * € L 2 
provided (3 2 < S — £7p T . Applying Proposition 13.21 to the constant minimizing sequence 
= for which the sequence of densities p k — p obviously is concentrated, we see that 

> E^ T + AE N {1 - J E )-0{y/e), (51) 

where J £ is compactly supported and e is small enough. Since e can be arbitrarily small 
we obtain that E > Ep T + AE N , which concludes the proof. □ 



5 Proof of Binding 

In this Section we prove Theorem 12.11 part (b) . To explain the main ideas in their pure 
form, without the difficulties due to the Pauli-principle, we first do the proof in the case 
of Bolzons, i.e., for Pekar-Tomasevich functional defined on L 2 (K 3W ). Thereafter we shall 
describe the modifications necessary to accommodate fermions and bosons. 

The case of Boltzons. The functionals £ x , . . . ,£ N ~ 1 have minimizers $i,...,$jv-i by 
assumption. Assuming that 

E prp = Eprp — Eprp (^2) 

for some k £ {1, ...,N — l}we shall prove in the Steps 1 and 2 below, that on the one hand 
&k ® &N-k is a minimizer of £ N , on the other hand it cannot satisfy the corresponding 
Euler-Lagrange equation. Hence the assumption (|52[) must be wrong. 

Step 1: $fe ® &N-k is a minimizer of £ N , that is 

£ N (*k®*N-k)=EpT' (53) 

From the definitions of the density, p$, and interaction energy D(p$) associated with 
any "J (see ©, ©), we easily see that 

N — k 

(54) 

and 

D(p<s> k ®<s> N _ k ) = D(p* k ) + D(p* N _ k ) + 2D(p* k , p$ N _ k ), (55) 

where 

k N ^ 

2D(p< s > k ,p< s , N _ k ) = (<f> k ®$ N _ k ,J2 E — r^k^^N-k)- (56) 

i=l.J=fc+l lafi ^'1 

From (|55|) . ([56]) . and the assumption (|52|) it follows that 

(8 $Ar_ fc ) = £ k (<S> k ) + £ N - k ($ N - k ) 

Ej prp -\- Ej prp Hipp. 

Step 2: $fc ® $jv-fc does not solve the Euler Lagrange equation of £ N . 
Suppose that <g> $jv-fc solves the Euler-Lagrange equation 

N N N \ 

E D U + E uTTFi " S W**-* - A J $fe ® = 0, (57) 

for some A G 1. Since $ k and 3>N-k are minimizers of £ k and £ N ~ k , respectively, they 
satisfy the Euler-Lagrange equations 

E + E T^-—n - E ^ - A i ** = o. ( 58 ) 
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and 

N N N , 

E D U + E - E Vp.„-> (%) - A 2 ] = °' ( 59 ) 

j=k+l k+l<i<j ' 1 Jl j=fe+l 



with Ai, A 2 G M. Note that, by (J5] 

Taking tensor products of the Equations (|58p and (15^1) with $jv-fc and 3>fc, respectively, 
and subtracting the resulting equations from (1571) . we obtain that 

fe N k n \ 

E E l3h"EV»W- E ^fe)-A + Ai + A 2 U fe ®$ JV _ fe = 0. 

i=l j'=fe+l 1 1 Jl i=l i=*+l / 

(61) 

Since V p & and V^„, n fc are bounded functions (see (|104[) in the Appendix) the expression 
in parentheses is a multiplication operator that is bounded below by 

k N 

E E -r~r M > (62) 

for some M > 0. Clearly, (|62|) is positive, e.g., for cci close to Xk+i- We may thus find balls 
Bi C K 3fe and B 2 C R 3 ^"^ such that {62]) is strictly positive on Bi x B 2 . At the same 
time we may assume, after suitable magnetic translations of <&N-k, that 

|$fc®*w-fc| 2 >0. (63) 

B, xB 2 



The strict positivity of the lower bound (|62[) and and the inequality (|63[) are in contradiction 
with ([BT]) . which completes the proof of Step 2. 



The case of fermions. In the case of fermions, the tensor product (g> $at_a; of the 
minimizers and •J'jv-fc in Step 1 must be antisymmetrized and normalized. The density 
of the resulting TV-particle state is not the sum of the densities of 3? at-/,-. In order to 
regain an analogue of (|54")1 we shall apply smooth space cut-offs at distance R from the 
origin and then move <5>N-k by a distance of SR. These cut-off minimizers, as well as 
their antisymmetrized tensor product, are approximate minimizers satisfying approximate 
Eulcr-Lagrangc equations, the error being exponentially small. But such an exponentially 
small error is not compatible with the power laws decay of the Coulomb interaction between 
the first k and the last N — k particles. 

We now proceed with the details. We use c, C to denote positive constants possibly 
changing from one equation to another. Suppose that for some k 

Ep T = Ep T + Ep T fe , (64) 

and let ip m be a minimizer off" 1 , m G {1, . . . , N- 1}. Let / G C°°{R; [0,1]) with f(s) = 1 
if s < —1 and f(s) = if s > 0, and let xr( x ) '■= f(\ x \ ~ R), a smoothed characteristic 
function of the ball B(0, R) C M 3 . We define 

, _ tor 

\WmXTW 

Let y G M 3 with |y| = 3R. Recall that T^jv-fc denotes a magnetic translation of <pN-k as 
defined in (lll[) . Due to the exponential decay of the minimizers ipk,^N-k and their gra- 
dients and Laplacians we obtain that <f>k,Ty<f>N—k are approximate minimizers and satisfy 
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respectively the Euler-Lagrange equations of £ k ,£ N k up to an exponentially small error. 
More precisely, 

{H k pH ,-E k PT )^ k = Qk , A {e- cR ), (65) 
( H PT~X- h ~ E PT k )Ty<t>N-k = QN _ k:A (e~ cR ), (66) 
where Og m A refers to the Q m ,A norm. Since £ m (<fr) — (cj),H™({>) it follows that 

£ k (M = E k T + 0(e- cR ), 
£ N - k (T y ct>N-k) = E»T k + 0(e-° R ). 



Equations (|65|) and (|66j) correspond to (|58[) and (f59|) in the boltzonic case, note however 
the irrelevant constants D{p r j )k ) and D(pT y( p N _ k ) in the Hamiltonians defined by (|2ip . 

Let now $ := P k (<p k ®T y 4>^^ k ). Here := \J~(^)Pa where Pa denotes the projection 
onto the completely antisymmetric functions with respect to permutations of pairs of posi- 
tions and spins. The factor in front of Pa is chosen so that $ is also normalized. Since the 
densities of <fi k ,T y (f>N^ k have disjoint supports we obtain that p$ = p^ h + pT y cf> N _ k which 
similarly to the case of Boltzons implies that 

S N {$) =£ k (M +£ N - k {T y $ N - k ) (68) 

and that 

V P *=V P4k +V PTy , N _ k . (69) 
From (|67|) and (|68|) we obtain that 

£ JV ($) = E k PT + E^ k + 0{e- cR ). (70) 

We show now that <I> satisfies an approximate Euler Lagrange equation. We take the tensor 
product of both sides of |65|) with T y cf)N-k- Similarly, we take tensor product of both sides 
of (|66|) with <fi k . By adding the resulting equations and then adding Jk{4>k <B> T y (f>N- k ) on 
both sides, where 

k N ^ k N 

j k ( Xl ,...,x N )-.= j2 \ x .- x .\ ~Ew t ^)" v P*„(. x i)> ( 71 ) 

i=l j=k+l 1 1 3 ' i=l j=k+l 

we arrive at 

~ c R)<t>k ® T y( j> N _ k = J k ^k ® P,0iv-fc) + O s „ iA (e- cit ), (72) 

where := Pp T + Ep^ k — 2D(p t j >k , PT v <p N _ k ) depends on i?. We have used (j2"T1) and 
(|69p . The fact that the supports of 4>k, T y 4>N-k have distance R in each particle coordinate 
implies that 

\J k \ = 0(R- 1 ), \VJ k \=0{R- 2 ), uniformly in (x u ...,x N ) & supp0 fc ®T y <j> N - k . (73) 

Applying the antisymmetrization Pk to both sides of (|T2[) and using (1751 as well as the 
symmetry of Hj^ with respect to the N particles, we arrive at 

KHZ-CRMa^^OiR- 1 ). (74) 

We are now going to improve this error estimate by changing the Lagrange multiplier by 
0(R~ 1 ). To this end we write 

$ = A fl $ + f R , with (f R , $) = 0. (75) 
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First observe that (|74|) and ff75|) imply that Xr = cr + 0(R 1 ) and therefore 

ll^lls^^OCi?- 1 ). (76) 

On the other hand using ([75)1 twice we obtain that 

(/«, /«) = (fa, H»*) = (f R , « - E» T )$). 

Recall that (Hj^ - Ep T )\ UN > (see Lemma l2^4|) and e Hjv- Hence, by Cauchy- 

Schwarz for positive (semi-)definite quadratic forms we find 

ll/flll 2 < </«, « ~ ^ r )/H> 1/2 <$, « - E$ T W* 

<4fR\\Q»A^iK- E PTW 2 - 

This estimate together with = ($,ff^$) (see Lemma l2l|. (|64| . ((7Qj) and ([76]) 

implies that 

II(^-Ak)$|| = ||/ii||=0(e- oii ). (77) 
By definition of $, the equations (jT5l) and ([77| imply that 

|| (^ - Afl^fcgT^fcll - 0(e- cR ), (78) 

because acts isometrically on the left hand side of ([78")) and commutes with . From 
(|72|) and <j78j) it follows that 

||(J fc + c fl - A R )<fe ® T„^_ fc || = 0(e- cif ). (79) 

This is in contradiction with Lemma l5.1l below. Hence, our assumption (164[) must be wrong 
and (fTUf is proved. 

Lemma 5.1. If the minimizers tp m , m € {fc, TV — fc} are chosen so that J up^ m {u)du = 0, 
then there exists a constant C > smc/i that 

C 

uif R ||(J fc + M)0 fe ® Ty^.feU > -^3 
(recall that \y\ = 3i?j. In particular, (|79j) does noi ZioZd. 

Proof. Let M £ Rbe arbitrary. Recall that Tj,</>jv_fc by definition is a magnetic translation 
by y with |y| = 3i? of 4>N-k- By a change of variables for the particles with labels in 
{k + 1, ...,N} we find that 

IKJfc + AO&ST^w-fcll =||/h||, (80) 

where 

/if := (Jk + M)<j) k <g) ^jv^fe (81) 

and 

k N ^ k N 

J k ( Zl ,...,z N ):=J2 E u._,. +1/| -£W*^ + y)- E Vp*(*i-V)' ( 82 ) 

i=l j=fc+l 1 1 3 1 »=1 J=fc+1 

By ([80]) it remains to prove that there exists C > independent of M so that 

\\Ir\\ > |p (83) 
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From the assumption of the lemma and the exponential decay of ip m we obtain for <j} m that 
u P<t , m {u)du = 0(e~ cR ), for m G {fc, N — k}. (84) 
By normalization of ^ m and the definition of 4> m we may choose d > such that 

/ dz 1 ...dz N \(t) k (z 1 ,...,z k )\ 2 \(t) N - k (z k+1 ,...,z N )\ 2 > i (85) 

for all i? > d + 1 where B d := B(0,d) N C R 37V . To prove <j83j) , and thus the lemma, it 
clearly suffices to show that 

\\lRXB d \\>§z, C>0, (86) 

where C is independent of M. 

We are going to expand ([82]) in powers of t^t . To this end we first remark that 

i , m-™? - m 2 , ( Hiy uniformly in H < 2i?; (87) 



1^ — 2/1 |y| \y\ 2 2|y|3 VN 4 

where y = y/\y\- Using this, (|54"|). suppp^ m C B(0,R) and the definition of U p (see (|2U)) ) 
we obtain for \zi\, \zj\ < d 

Vf ,_ ( « +9) = <* - *, - ^ + ^£fM! ) + + 0(fl - 4) , (89) 

where 

JmW := / A/> m (") ^ (90) 

Recall that (j) m depends on R and hence the exponential decay of ip rn is needed for estab- 
lishing the bound 0(i?" 4 ). Using flST]) again we obtain 

1 (1 y-(zj-Zi) 3(y ■ (zj - Zi)) 2 - \zj - z.^ 2 



\zi-Zj+y\ \\y\ \y\ 2 2\y\ 3 

+ 0(R- 4 ), V^j, Zj G B(Q, d). (91) 

Inserting (f88j) . (l89l) and (l9Tj) into (IBTl) (see also ([82]) ) an elementary but somewhat lengthy 
calculation gives that 

1 I k N \ 

I R XB d {zi, ■ ■ ■ , z N ) = — 2J ^ ( 2i • Zj - 3(z; • y)(zj • y)) + C y . M 
m \i=ij=fe+i / 

x (0 fc ® </> N -k)XB d ( z i> -,z N ) + 0(R- 4 ), (92) 

where C„,m = -(iV-fc)/fe(t/)-fc/Ar_fc(y)-/c(iV-fc)|y| 2 + Af|y| 3 depends on y and M only. 
We recognize in (|9"2"|) the interaction energy (z^ • zj — 3(zi ■ y)(zj ■ y))/\y\ 3 of two dipoles Zi 
and Zj separated by y. Let 



k N 



L (y> D ) = I dz x dz 2 ...dz N ^2 ( Zt ' z i ~ 3 ( z ' ' y ^ z i 'V))+ D 

jBd i=l j=k+l 

x (f>k(zi, .-, z k )(j) N -k{zk+i> •••) In) 



(93) 
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Then, by ([92]), 

1 

w 

and it remains to show that there exists a constant C > such that 



IrXbA = T^Hy,C y , M ) 1/2 + 0(R- 4 ), (94) 



L(y,D)>C, Vy,D. (95) 



This estimate together with (|94l) concludes the proof of (|86| and therefore of Lemma 15.11 
To prove (j9"5")) we first establish that L(y, D) is everywhere positive. To this end we fix 
y, D and we consider the function 

f(zi,z 2 ,...,z N ) = ( X!*) ' ( X z j J ~ 3 ( ( X ^-yJ+A 

which is part of the integrand in (|9"5)l . One can show that /(zi, Z2, zjv) ^ almost 
everywhere, which together with (|85p implies that 

L(y,D)>0. (96) 



Now we will use a continuity argument to show (|95[) . Since / is continuous and thus 
bounded on Bd it follows, by the dominated convergence theorem, that L is a continuous 
function of y, D. Moreover, Ummi^oo L(y, D) = oo. These observations together with (|M|) 
and the fact that a continuous function on a compact set attains its minimum give ([55)1 . □ 

A Properties of E N 

Lemma A.l. 

(a) The functional £ N is bounded from below on the set Sm '■= S Qn,a '■ ll^ll = !}■ 

(b) Every minimizing sequence (^k) of £ N on Sn is bounded in Qn,a- 

(c) If for a minimizing sequence (^k) of £ N we have that — > ^ weakly in Qn,a and 
strongly in Hn then "J is a minimizer of £ N . 

Proof. (a),(b) They follow from the Hardy and diamagnetic inequalities. We recall from 
[T4"] that the diamagnetic inequality states that if <j> € H\ then we have that |V|</)|(x)| < 
\(D A, Xl <j){x) , Da, xn <P(x))\, for almost all x £ K 3Ar . In the case N = 1 and without spin 
a detailed proof of parts (a) and (b) of the Lemma is given in [8] and in the general case 
the argument is similar. 

(c) We will now show that if "ffc — > *f> weakly in Qn,a and strongly itiHn then 

£ N (V) < liminff^fe), (97) 

k— >oo 

from which we conclude that ^> is a minimizer of the Pekar-Tomasevich functional £ N . 
Indeed, recall that 



e»m = {% (f; D% Xj + x -±- ) *> - do*). (98) 



Since * weakly in Qn.a and \\^k\\ = 11*11 = 1 we see that 

N N 



(*,E£»^.*) <liminf(* fc ,E£»i )X .* fc ). (99) 

3=1 J'=l 
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On the other hand, since ^> k — > ^ weakly in Qn,a and since \xi — Xj\ 1 is a bounded 
operator from Q n,a to Hn we obtain that 

\xi — Xj \~ ll i'k — > \xi — Xj\ ^f, weakly in Hn- 

Since, moreover, ty k — > ^ strongly in Hn we conclude that 

\x z - x.r 1 *) = lim (* fc , \ Xi - 1,-r 1 **). (100) 

k— >OQ 

In addition, 

D(jy<st) - D(p k ) = D(p^ - p k ,pq,) + D(p kl p^ - p k ). (101) 

We will show that 

D(p k ,p*-p k )^0. (102) 
Indeed, using (|SJ) and (j2"0")) we obtain that 



D(p k , py - p k ) = J V Pk (py - p k )dx. (103) 

But using Lemma [A.. 11 (b) together with (|20p. (|3"T|) with r = oo, and (j3"2")l we can prove that 

sup||VpJ| L - < oo. (104) 

k 

Since ^ k —> 'J in Hn we obtain that — p^-Hl 1 which together with (jl03[) and 
(fT04| implies (fj"02j) . Similarly, 

-D(p*-Pfc,P*)^0. (105) 
Combining lfT02]> . (fT05|) and (flUTjl we obtain that 

£>(p*) = lim D( Pk ). (106) 
The relations ([98 ]) . (p f , (fTOO)) and (fT06|) give (|97]) as desired. □ 

References 

[1] Shmuel Agmon. Lectures on exponential decay of solutions of second-order elliptic 
equations: bounds on eigenf unctions of N-body Schrddinger operators, volume 29 of 
Mathematical Notes. Princeton University Press, Princeton, NJ, 1982. 

[2] Ioannis Anapolitanos and Benjamin Landon. The ground state energy of the multi- 
polaron in the strong coupling limit. arXiv.1212.3511 . 

[3] F. Brosens and J. T. Devreese. Stability of bipolarons in the presence of a magnetic 
field. Phys. Rev. B, 54:9792-9808, Oct 1996. 

[4] F. Brosens, S. N. Klimin, and J. T. Devreese. Variational path-integral treatment of 
a translation invariant many-polaron system. Phys. Rev. B, 71:214301, Jun 2005. 

[5] H. L. Cycon, R. G. Froese, W. Kirsch, and B. Simon. Schrddinger operators with 
application to quantum mechanics and global geometry. Texts and Monographs in 
Physics. Springer- Verlag, Berlin, study edition, 1987. 

[6] Rupert L. Frank, Elliott H. Lieb, Robert Seiringer, and Lawrence E. Thomas. Stability 
and absence of binding for multi-polaron systems. Publ. Math. Inst. Hautes Etudes 
Sci, (113):39-67, 2011. 



18 



[7] G. Friesecke. The multiconfiguration equations for atoms and molecules: charge quan- 
tization and existence of solutions. Arch. Ration. Mech. Anal., 169(1):35-71, 2003. 

[8] M. Griesemer, F. Hantsch, and D. Wellig. On the magnetic Pekar functional and the 
existence of bipolarons. Rev. Math. Phys., 24(6):1250014, 13, 2012. 

[9] M. Griesemer and D. Wellig. Strong coupling polarons in electromagnetic fields. In 
preparation. 

[10] Marcel Griesemer and Jacob Schach M0ller. Bounds on the minimal energy of trans- 
lation invariant AT-polaron systems. Comm. Math. Phys., 297(l):283-297, 2010. 

[11] M. Lewin. Geometric methods for nonlinear many-body quantum systems. J. Fund. 
Anal, 260(12):3535-3595, 2011. 

[12] Elliott H. Lieb. Existence and uniqueness of the minimizing solution of Choquard's 
nonlinear equation. Studies in Appl. Math., 57(2):93-105, 1976/77. 

[13] Elliott H. Lieb. Bound on the maximum negative ionization of atoms and molecules. 
Phys. Rev. A, 29:3018-3028, Jun 1984. 

[14] Elliott H. Lieb and Michael Loss. Analysis, volume 14 of Graduate Studies in Mathe- 
matics. American Mathematical Society, Providence, RI, second edition, 2001. 

[15] Elliott H. Lieb and Walter E. Thirring. Universal nature of van der waals forces for 
coulomb systems. Phys. Rev. A, 34:40-46, Jul 1986. 

[16] P.-L. Lions. The concentration-compactness principle in the calculus of variations. 
The locally compact case. I. Ann. Inst. H. Poincare Anal. Non Lineaire, 1(2):109- 
145, 1984. 

[17] G. Verbist, F.M. Peeters, and J. T. Devreese. The stability region of bipolaron forma- 
tion in two and three dimensions in relation to high-tc superconductivity. Solid State 
Communications, 76:1005-1007, Nov 1990. 



19 



